В Тимано-Печорской нефтегазоносной провинции установлен широкий стратиграфический диапазон залежей нефти и газа, приуроченных к отложениям ордовика, силура, девона, карбона, перми и нижнего триаса. Тектонической этажностью обособляются формационные карбонатные – ордовикско-нижнедевонский, верхнефранско-турнейский, верхневизейско-нижнепермский и терригенные – среднедевонско-нижнефранский, нижневизейский и уфимско-триасовый нефтегазоносные комплексы. Ни в одном районе СНГ нет такого многообразия спектра нефтяных, газоконденсатных и газовых скоплений, контролируемых разными формами ловушек, несущих черты комбинированного сочетания структурного и неструктурного генезиса. В этих условиях рациональное комплексирование геологических и геофизических методик при поисках и разведке нефтяных и газовых месторождений является приоритетной задачей научного исследования.

Цель выполненных нами исследований – определение емкости месторождений в межскважинном пространстве по данным сейсморазведки (СЕМП).

Исследования по данной проблеме осуществлялись на двух объектах Хорейверской впадины, размещенных в глубокозалегающем и наименее изученном и, как до недавнего времени предполагалось, малоперспективном ордовикско-нижнедевонском комплексе осадочного чехла. Один из них – месторождение им. Р. Требса (Варкнавтская площадь) – приурочен к отложениям нижнего девона северо-восточного борта, другой – Верхневозейское месторождение – к отложениям нижнего силура юго-западного борта. Вмещающие образования – доломиты и известняки, региональная покрышка – глинистый кыновско-саргаевский экран. Оба месторождения крупные, контролируются комбинированными структурно-литологическими ловушками.

Месторождение им. Р. Требса по продуктивным отложениям не имеет ясных очертаний. В рельефе нижнего пласта Д₁¹ изогипсой минус 4100 м с запада оконтурен 7-километровый полувал меридионального простирания с четким западным склоном, очерчивающим два локальных купола амплитудой по указанному направлению около 50 м. Их восточное ограничение не замкнуто.

Верхний пласт Д₁² залегает в ложбине, развитой на месте полувала с тальвегом – 4000–4100 м и бортовыми уступами высотой до 120 м. Нефтяные залежи приурочены, в основном, к карбонатным толщам лохковского яруса нижнего девона, представленным разрезами двух типов. В скв. 1, 5, 7 это доломиты и известняки мелкокристаллические, крепкие, трещиноватые, местами нефтенасыщенные, фиксируемые минимальными значениями ГК (до 28·10^-14 А/кг) и повышенными HFK (5–7 усл. ед.). Восточнее, в скв. 6, 9, 12, породы-коллекторы замещены уплотненными глинистыми карбонатами с повышенными значениями ГК, ИК и пониженными НГК. Региональной покрышкой для залежей нижнего девона служат кыновско-саргаевские отложения франского яруса мощностью 83–156 м.

Верхневозейское месторождение приурочено к одноименной структуре, прилегающей к Колвинскому авлакогену. Размеры структуры по подошве кыновских отложений составляют 28x11 км, амплитуда – до 250 м в пределах изогипсы минус 3480 м. В юго-восточной части к основной структуре примыкает небольшое куполовидное поднятие (район скв. 216) с амплитудой около 50 м и размерами 4,6x2,5 км. Севернее выделяется небольшое тектонически ограниченное поднятие (район скв. 218) с амплитудой 30 м и размерами 0,8x4,8 км. Отмечаются детальные сбросовые нарушения на сейсмопрофилях в районе скв. 205, 206, 207, 234, 201, 211 со смещением блоков до 60 м.

Залежи нефти на Верхневозейском месторождении установлены в доманиковых и нижнесилурийских отложениях. Залежи в отложениях доманикового горизонта имеют линзовидный тип, характеризуются АВПД и залегают на глубинах 3279–3557 м. Залежи нефти в отложениях нижнего силура приурочены к отложениям верхнемакарихинской толщи сандивейской и веякской свит.

Основная залежь в нижнем силуре Верхневозейской структуры приурочена ко всем выделенным свитам, является пластово-массивной, тектонически экранированной. Размеры залежи 16x4,5 км, этаж нефтеносности около 260 м с ВНК. на отметке минус 3484 м (по скв. 210, 204, 207 м). Глубина залежи 3290–3550 м.

В отличие от традиционной методики прогнозирования геологического разреза (ПГР), использована методика количественного определения параметра эффективной удельной емкости коллектора q в межскважинном пространстве по данным сейсморазведки МОГТ [2]. Эта величина, как известно, представляет собой произведение средневзвешенного значения коэффициента пористости коллектора К_п и его суммарной толщины Sh_эф, установленных по данным глубокого бурения.

Результаты сейсморазведки МОГТ, обработанные на ЭВМ, дают возможность заменить обычно интерполируемые значения параметра q экспериментальными значениями по данным обработки фактических сейсмических измерений.

Привязка этих данных к геологическим разрезам, ГИС, сейсмокаротажу (СК) и акустическому каротажу (АК) должна показать их связь с потенциально продуктивными отложениями. Значения q в межскважинном пространстве, в конечном счете, определяются по изменениям временного интервала низко- и высокоскоростной частей ПАК (Dt_пак) или псевдоакустических скоростей (DV_пак).

Применяемая методика определения параметра в межскважинном пространстве базируется на результатах отечественных и зарубежных исследований [1, 3–5] детальных связей между акустическими свойствами и петрографическим характером карбонатных пород в широком диапазоне фаций, пористости, литологического состава, вида скелета породы и показавших, что главным фактором, определяющим импеданс и скорости продольных и поперечных волн, является пористость. Все прочие переменные, включая первичный состав, вторичные литологические изменения, петрофизический состав скелета, тип порового флюида не оказывают заметного влияния на величину скорости. При этом продуктивные карбонатные отложения характеризуются пониженными значениями акустических и псевдоакустических скоростей, а, следовательно, повышенными значениями средних коэффициентов пористости.

В процессе определения q используются не псевдоакустические скорости в коллекторе, зависящие, главным образом, от Кп, а Dt_пак и DV_пак, соответствующие продуктивным отложениям, содержащим коллекторы и неколлекторы и функционально связанные как со средней пористостью коллектора (K_п), так и их долей в общей толще продуктивных отложений (долей суммарной эффективной толщины коллекторов Sh_Эф).

Таким образом, при увеличении Кп скорость уменьшается, а Dt увеличивается. То же происходит и при неизменном К_п, но при увеличении общей толщины коллекторов Sh_эф. Когда исследуется DV_пак B ЭКО-пласте, соответствующем продуктивным отложениям, при увеличении Кп V_пак в пласте уменьшается, a DV_пак, т. е. относительные изменения V_пак увеличиваются. То же происходит и при увеличении Sh_эф, но здесь все зависит от выбора уровня, относительно которого определяются DV_пак: если DV_пак рассчитывается относительно минимальной скорости, соответствующей максимальной пористости, то увеличение DV_пак происходит с уменьшением Кп.

При изучении массивных залежей УВ используется Dt_пак, а для характеристики параметра удельной емкости коллектора залежей пластового типа –DV_пак.

Функциональные зависимости Dt_пак= (K_пSh_эф) или DV_пак= (K_пSh_эф) определяются экспериментально по данным бурения, ГИС и сейсморазведки (ПАК, ЭКО) для конкретных месторождений или групп однотипных месторождений.

В общем виде эти зависимости обычно представляют собой полином невысокой степени: Dt_пак, DV_пак= Aq, где А – оператор корреляционной связи. Значения Dt_пак, DV_пак устанавливаются по разрезам ПАК и ЭКО, и по зависимости Dt_пак, DV_пак=f(K_пSh_эф) пересчитываются на сейсморазведочных профилях в межскважинном пространстве в значения q. Таким образом, интерполяционные значения q заменяются в межскважинном пространстве расчетными по данным сейсморазведки.

Для построения емкостных моделей продуктивных карбонатных нижнедевонских и нижнесилурийских отложений соответственно на Варкнавтском и Верхневозейском месторождениях проведена специальная обработка сейсморазведочных данных МОГТ на ЭВМ ЕС и ПС-2000, заключающаяся в получении временных разрезов МОГТ с восстановлением (сохранением) истинных амплитуд и последующей их трансформации в разрезы ПАК, ЭКО и частично в псевдоакустической жесткости (ПАЖ) Исходя из задач количественной интерпретации, во временных интервалах регистрации целевых нижнедевонских и нижнесилурийских отражений была сделана подискретная распечатка значений псевдоакустических скоростей по трассам ПАК для каждого из профилей с шагом t₀ = 4,8 мс и х= 100–500 м.

Вся обработка на ЭВМ приведена с использованием отечественных систем СЦС-3-ПГР (ЦГЭ МНП), СОС ПС и ПАРМ (ВНИИгеофизика), РЕАПАК (СНИИГиМС) на представительном объеме – 40 профилей длиной 644 км (Варкнавтская площадь 17 профилей, 309 км, Верхневозейская 23 профиля, 335 км).

Кроме обработки сейсморазведочных профилей МОГТ проведена обработка данных акустического каротажа по 20 скважинам (7 скважин на Варкнавтской площади и 13 скважин на Верхневозейской) с построением графиков зависимостей пластовых акустических скоростей (V_пл^АК) от глубины (H), а в тех скважинах, где проведен сейсмокаротаж (СК) – зависимостей V_пл^АК=f(t₀) для прямого сопоставления акустических и псевдоакустических параметров изучаемых отложений.

На Варкнавтской площади на ЭКО-разрезах выделены и прослежены практически непрерывно два сейсмических пласта, соответствующих высокоскоростным продуктивным нижнедевонским геологическим пластам Д₁² и Д₁¹.

Величина временной невязки данных СК, АК и ЭКО-, ПАЖ-разрезов не превышает 7–10 мс, что гарантирует высокую точность последующей интерпретации.

На большей части Верхневозейской площади, где выделены все четыре нижнесилурийских пласта, первые два сверху (1C и 2С) соответствуют веякско-сандивейским отложениям, а нижние (3С и 4С) – сандивейско-макарихинским и салюканским отложениям.

Сейсмические пласты 3С и 4С прослежены по всей площади, а 1C и 2С выклиниваются в общем направлении с востока на запад, что отождествляется с изменением акустических свойств нижнесилурийских отложений в связи с выклиниванием веякских пород вблизи Колвинского мегавала.

Разрезы ПАК дали возможность изучить общие закономерности изменения псевдоакустических скоростей в пластовой дифференциации.

Дискретные распечатки V_ПАК позволили охарактеризовать псевдоакустическими скоростями выделенные на ЭКО-разрезах сейсмические пласты

Для проведения дальнейшей интерпретации необходимо было убедиться, что характер изменения V_АК и V_ПАК в сейсмических пластах одинаков, поскольку при интерпретации ПАК-преобразова-ний использованы относительные изменения V_ПАК по латерали – DV_ПАК которые должны изменяться так же, как и DV_АК.

На рис. 1 приведен результат примерно адекватных изменений DV_ПАК и DV_АК в районе глубоких скважин на Верхневозейской площади. Подобные соотношения служат основанием для использования значений DV_ПАК в межскважинном пространстве и их пересчета в параметр удельной емкости коллектора q=K_пSh_эф по выявленным статистическим зависимостям между этими параметрами для каждого пласта.

На рис. 2 показана установленная статистическая связь между DV_ПАК и параметром q=K_пSh_эф для одного из пластов. Подобные графики построены для всех исследуемых пластов. Эти зависимости были использованы для пересчета DV_ПАК в q = K_пSh_эф в межскважинном пространстве, т. е. трансформации карт DV_ПАК в карты q.

Рассмотрение новой геолого-геофизической информации об емкостных свойствах продуктивных отложений нижнего девона (рис. 3) и силура (рис. 4) свидетельствует, прежде всего, о наличии на обеих площадях существенно нелинейного закона изменения параметра удельной емкости с образованием локальных максимумов и минимумов. Такое распределение емкостных свойств девонских и силурийских продуктивных отложений на Варкнавтской и Возейской площадях объясняет возникновение больших ошибок при любой межскважинной интерполяции параметра q = K_пSh_эф или просто К_п.

1. Методология и некоторые результаты комплексного изучения сложно построенных геологических сред методами ГИС и сейсморазведки при поисках и разведке месторождений нефти и газа / О.Л. Кузнецов и др. // Новые геоакустические методы поисков и разведки месторождений полезных ископаемых // Тр. ВНИИЯГТ –М, 1986
2. Определение параметра удельной емкости коллектора в межскважинном пространстве / Е.А. Копилевич и др. // Геология нефти и газа.– 1988– № 8. С 37–41.
3. Dominica S.N. Rock lithology and porosity determination from shear and compressional wave velosity – Geophysics, vol 49, N 8, 1984
4. Got H., Mirable I. Seismic hthology is one methods of sedimentology – Geophysics, vol. 49, N 10, 1986.
5. Ratavtch F., Kendall С.H. , Todd Т.F. The relationship between acoustic properties and the petrographic character of carbonate rocks – Geophysics, vol. 49, N. 10, 1984.

Effectiveness of the use of the method of quantitative determination of resorvoir holding capacity in inter-well space by MCDP seismic data is shown at the example of the Devonian and Silurian carbonaceous deposits study at Varknavtskaya and Verkhne-Vozeyskaya areas of Timan-Pechora province. Possibility of DV_PAL parameter use (lateral increasing in beds) instead of usually used V_PAL parameter (interval beds velocity of pseudoacustic model) is regarded for these purposes.

Рис. 1. Изменения DV_пак и DV_ак для пласта 3с в районе скважин на примере Верхневозейской площади

Рис. 2. Зависимость DV_ак и DV_пак от q=K_nSh_эф для пласта 3с на примере Верхневозейской площади:

1 – коэффициент пористости по данным АК (Кп^АК ); 2 – коэффициент пористости по данным НГК (Кп^НГК ); 3 – дисперсия К_пh_эф; 4 – номера скважин

1 – изолинии равных значений параметра удельной емкости коллектора; 2 – предполагаемые контуры запасов залежей; 3 – рекомендуемые и пробуренные скважины; 4 – глубинный разлом по данным сейсморазведки

1 – изолинии равных значений параметра удельной емкости коллектора; 2 – контуры запасов залежей; 3 – скважины